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Abstract The quantum behavior of a precooled cantilever can be probed highly efficiently 

by electrostatically coupling to a trapped ultracold ion, in which a fast cooling of the cantilever 
down to the ground vibrational state is possible. Within a simple model with an ultracold ion 
coupled to a cantilever with only few vibrational quanta, we solve the dynamics of the coupling 
system by a squeezed-state expansion technique, and can in principle obtain the exact solution of 
the time -evolution of the coupling system in the absence of the rotating-wave approximation. 
Comparing to the treatment under the rotating-wave approximation, we present a more accurate 
description of the quantum behavior of the cantilever. 

PACS numbers: 42.50. Vk, 61.46.+W, 62.25. +g 



I Introduction 

With recent rapid progress in microfabrication technology, the quantum electromechanical system 
(QEMS), with nano- to micro-meter scale mechanical resonator electrostatically coupled to 
comparably sized electrical devices, has been drawn much attention [1]. Typical QEMSs include 
displacement detectors strongly coupled to high-quality cantilevers or doubly-clamped beams 
[2,3], or to single electron transistors [4,5]. With current state -of-the -art techniques, the quantum 
behavior of the QEMS can also be efficiently evidenced by the coupled electrical devices under 
the proper temperature and vacuum conditions. For example, detection of quantum superposition 
of position states of mechanical resonators has been available by using an electrostatically coupled 
Cooper-pair box [6]. In addition, it is expected that single phonon detection [7] and quantum 
tunneling of mechanical degrees of freedom [8] could be done by means of QEMS. Moreover, 
with the idea of QEMS, it is becoming an exciting research frontier to couple atomic systems to 
nano-mechanical devices for demonstrating some new phenomena in condensed matter physics 
[1,9]. 

Besides its own interest, QEMS is also playing more and more important role in the context of 
quantum information processing (QIP). We have noticed an idea [10] for QIP with coupled 
trapped ion to a high-quality nano-mechanical resonator. Under certain conditions, the two 
component subsystems interact and yield energy exchange, entanglement and quantum gating. 
This is a promising way towards large-scale QIP. Moreover, readout of individual qubits is very 
challenging in solid-state QIP, which is resulted from the difficulty of the single spin detection. A 
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recent breakthrough in this respect is the detection of a single unpaired electron-spin on the 
surface of silicon dioxide by magnetic resonance force microscopy [11], which is actually a 
QEMS. Another QEMS - single electron transistor based on a fullerene [12] also provides a 
promising way toward readout of qubits in fullerene-based QIP [13]. 

In this paper, we focus on the dynamics of a coupling system involving a very small doubly 
clamped cantilever (- a nonomechanical resonator as mentioned in [10]) and a single trapped 
ultracold ion, which has originally been studied in [14]. Our work can also be applied to [10] since 
there are many things overlapping between [10] and [14]. One of the purposes in [14] is to highly 
efficiently probe the vibration of a cantilever precooled to 4 Kelvin, through an electrostatically 
coupled trapped ion initially cooled to the vibrational ground state, as shown in Fig. 1. Since the 
deviations of the cantilever and the ion from their equilibrium positions in the oscillation are 
comparable to the equilibrium separation between them, a second-order perturbative expansion of 
the Coulomb interaction term could yield the coupling of the position operators of the ion and the 
cantilever. Experimentally this can be done simply by applying an external bias voltage at an 
electrode of the cantilever to couple the motion of the cantilever to the vibrational state of the ion. 
Since the detection of the quantum vibration of the trapped ion is a sophisticated work by optical 
methods - radiation with blue or red detuned laser [15] followed by electric shelving amplification 
[16], if the two component subsystems can be coupled and decoupled at will by adjusting the bias 
gate voltage, the scheme in [14] is available with current technique by following three steps: The 
first step is to switch on the coupling of the cantilever to the ion. At a suitable time point when the 
vibration of the cantilever is fully transferred to that of the ion, the coupling is switched off and 
the red or blue detuned laser will be on to couple the vibrational and the internal degrees of 
freedom of the ion. The last step is to read out the vibrational quanta by detecting the internal 
levels of the ion. 

The other purpose in [14] is to fast cool the cantilever down to the ground vibrational state. The 
numerical calculation in [14] has shown the highly efficient transfer of the phonon from the 
cantilever to the ion. So the cantilever reaches its ground vibrational state only after a half 
coupling period. 

The focus of the present paper is to exactly solve, within the simple model given in [14], the 
dynamics of the coupling system consisting of a cantilever and an ion both under decay due to the 
thermal environment, by a squeezed-state expansion technique. In Section II, we will present the 
analytical solutions by our technique for the systems under the rotating-wave approximation 
(RWA) and beyond the RWA, respectively. The numerical results demonstrated in Section III show 
the big differences of the RWA solution and the non-RWA one in the case of relatively strong 
coupling, which implies our investigation in the absence of the RWA to be essential to the 
detection and to the cooling of the cantilever. Some discussions about the validity and the 
experimental feasibility of our theoretical work as well as the conclusion appear in Section IV, and 
we put some detailed deduction steps in Appendix parts. 

II The solution by squeezed-state expansion technique 
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Starting from [14], we consider the first step by switching on the coupling between the cantilever 
and the trapped ion, which yields the following Hamiltonian in units of h — 1 , 

H = (G)-iT a )a + a + {v-iT b )b + b-K(a + a + )(b + b + ), (1) 

where go and v are vibrational frequencies of the cantilever and the trapped ion, respectively, and k 
is the coupling constant between the two vibational degrees of freedom. To make a general 
treatment, we remain the counter-rotating terms in Eq. (1). So our study from Eq. (1) could be 
made in both weak and strong coupling cases. T d and T b are the decay coefficients regarding the 
vibrations of the cantilever and the trapped ion, respectively. We assume a trial solution of the 
wavefunction 

\t) = p^^w^.Mo.v^ww | ^ = p{t) | A{t)) ^ (2) 

where p (t), a t (t), a 2 (t), and a 3 (t) are variables determined later. With Eq. (2), the Schrodinger 
equation of Eq. (1) yields, 

Ap(i) + p(i)d l (i)a + a + + p(t)d 2 (t)a + b + + p(t)d 3 (t)b + b + ]\ A.t)) 
= p(t){2(6)- iT a )a l (t)a + a + +((0- iY a )a 2 (t)a + b + 

+ (v-iT b )a 2 (t)a + b + + 2(v-iT b )a 3 (t)b + b + (3) 
- x[a 2 (t) + [la^a^f) + a 2 (t)]a + a + + [2a 2 (t)a 3 (t) + a 2 (t)]b + b + 
+ [1 + a 2 (t) + 4a 1 (r)a 3 (r) + 2a,{t) + 2a 3 (t)]a + b + ] }\ At)). 

Comparing the terms | A^t)), a + a + \ A^t)), a + b + \ A^t)), b + b + \ A^t)) in both sides of the equation, 



we obtain a set of equations, 

ip(t) = -Hra 2 (t)p(t), (4) 
id^t) = 2(o)-iT a )a l (t)-Ka 2 (t)-2KCc 2 {t)a l (t), (5) 
ia 2 (t) = (cq+ v-iY a -iT b )a 2 (t) - k- jca 2 (t) 

, Kb) 

- AKa x {t)a 3 {t) - 2ica x (t) - 2ra 3 (r)] 
id 3 (t) = 2(v -iT b )a 3 (t) - Ka 2 (t) - 2fca 2 (t)a 3 (t) . (7) 



From the initial condition CC 2 (0) = 0, (X 3 (0) = , and (t = 1 1 = 0) = 1, we have 

y(l-4|a,(0)| ! ) 

where the proof of Eq. (8) is put in Appendix I. Another initial condition is for the initial 
vibrational population of the cantilever, i.e., (t = 1«2 + «2| t = 0) = N a . With similar algebra to Eq. 
(8), we obtain 
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" (0)J -1-V . (9) 



i ,9 3/ a 

(1-4^(0) \y 2 



-V I N 

Eqs. (8) and (9) result in p(Q) = ±(N a + 1) /4 and a,(0) = ±- J - — ^ , from which we 

can straightforwardly obtain the exact dynamics of the coupling system by solving Eqs. (4)- (7). 
Before doing the numerics, however, for a comparison, we would like to solve the Hamiltonian in 
Eq. (1) under the RWA by the same technique we used above. For the Hamiltonian 

H = (0) - iT a )a + a + (v - iT b )b + b - K{ab + +a + b), ( 10 ) 

where the counter-rotating-wave terms are removed due to the approximation that rapid oscillating 
terms resulted by large detuning between the two vibrational degrees of freedom can be 
effectively averaged out in the calculation in the case of weak coupling. If we still assume Eq. (2) 
to be the trial solution of the wavefunction, we have the Schrodinger equation 

lp(t) + p(t)d l (t)a + a + + p(t)d 2 (t)a + b + + p(t)d 3 (t)b + b + ]\ A,t)) 

= p(t)[2(G) - iT a )a l (t)a + a + +{(0- iT a )a 2 (t)a + b + 

+ (v-iT b )a 2 (t)a + b + + 2(v-iT b )a 3 (t)b + b + (11) 
-2Ka x (t)a + b + -Ka 2 (t)b + b + 
-2fca 3 (t)a + b + - fca 2 (t)a + a + ]\A,t)}, 
which yields 

ip(t) = - p(t) = p(0) , (12) 
id x (t) = 2{w- jTJa.it) - jca 2 (t) , (13) 
id 2 (t) = (o)+ v-iT a -iT b )a 2 (t) - 2Ka x (t) - 2m 3 (f) , (14) 
id 3 (t) = 2(v -iT b )a 3 (t) - ica 2 (t) . (15) 

Different from the case in the absence of the RWA, we can find analytical solutions from Eqs. (13) 
to (15). After the direct algebra shown in Appendix II, we have 

a x {t) = a^e^' + < ) e £2<2V + <e n<3>f , 

a 2 (t) = a*y m ' + a i 2 l ) e aWt + a^' , (16) 
a 3 (t) = a^e am ' + a^e^' + <e a<3)f , 

where 

(1) _ _]_ I N a (2(Q-2iT a -iQ. (2) )(2Q)-2iT a -iQ. & ) 
a '° " (Q (3) -C2 (1) )(C2 (2) -C2 (1) ) ' 
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a (2) 



1 I N a (2v - 2iY b - iQ. (2) )(2w - 2iY a - ill (3) ) 



'" " 2\n 3 +i (Q (3) -a (2) )(a (2) -a (1) ) 



a 0) 



1 I N u {2(0 - 2iY a -iH {2) ){2v - 2iY b - ill (3) ) 
2\N a +1 (a (3) -a (1) )(a (3) -a (2) ) 



1 I iV a (2v-2/r 6 -ill m )(2(0-2iY a - ill (2) )(2(0-2iY a -£1 (3) ) 

- ! ° ~~2^sT7T (a (3) -a (1) xa (2) -a (1) )r " ' 

i I iv a (2v - 2/r 6 - n (2) )(2© - 2/r a - n (3) )(2<y - 2/r a - /ia (2) ) 



1*20 



2ViV a +l (Q (3) -C2 (2) )(ft (2) 



(17) 



a (3) 



1 J N a (2v - 2iY b - ill (3) )(2(0- 2iY a - ill (3) )(2eo- 2iY a - ill (2) ) 



20 2 AU+1 (q (3) -q (2) )(q (3) -q (1 V 



30 2\N a +\ (H (3) -H m )(H (2) -H m ) 



a (2) 



1 I N a (2(0 - 2iY a -iH {2) )(2(0 - 2iY a - ill (3) ) 



30 " 2AU +1 (Q (3) -C2 (2) )(Q (2) -Q (1) ) 



1 | N a (2(0-2iY a -i£l {2) )(2(0-2iY a -/1Q (3) ) 



a O) _ 

30 2 ^)iV a +l (C2 (3) -ft (1) )(£l (3) -C2 (2) ) 

III Numerical results 

We consider an achieved cantilever with vibrational frequency 19.7 MHz [3] and assume the 
frequency of the ion's vibration to be 16 ~ 19.7 MHz. For simplicity, we fix the decay coefficients 
T a , T b to be V 1 1000 . But the coupling strength K will be changed in our calculation to check 
the difference between the treatments under the RWA and by the non-RWA. Other initial 

conditions are n a (0) = N a = 6 and n b (0) = . In terms of our method presented in last section, 

no matter under which conditions, i.e., under the RWA or non-RWA, in order to obtain n a (t) and 

7i 6 (t) by solving 1 1) , we have to calculate (X x (t), OC 2 (t), (X 3 (t) and the initial conditions 
Eqs. (8) and (9). 

When the two vibrational degrees of freedom are nearly resonant, i.e., (0 ~ V , the difference 
between the treatments under the RWA and with non-RWA is almost invisible in the case of very 
small coupling constants (i.e., Klv < 0.1 in our calculation), as shown in Fig. 2. It can be well 
understood from the viewpoint of Jaynes and Cummings' seminal paper [17]: In this case, the 
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counter-rotating terms are negligible and the RWA can be safely used. Our result also confirms the 
validity of the RWA employed in [10] and [14] because Klv < 0.1 is satisfied there. Moreover, 
in this near resonant case, a weak coupling between the two vibrational degrees of freedom could 
achieve a quick cooling of the cantilever down to the ground state by a half coupling period, even 
in the case of relatively large decay coefficients T a and T b . 

To transfer the energy [10,14] or to make entanglement [10] with high fidelity, however, we prefer 
a large coupling K. While with the increase of the coupling strength (e.g, Klv > 0.1 in our 
calculation), the counter-rotating terms in the (near) resonant case is not negligible any more, and 
the RWA could not work well. This can be found in Fig. 3. To have a good transducer or to make 
high-quality quantum gates, we require the trapped ion to be decoupled from the cantilever at an 
exact time point. Therefore a non-RWA solution is essential to this situation. 

With the increase of/C, more and more differences between the RWA and non-RWA solutions 
would appear. Particularly, if the two vibrational degrees of freedom are detuned, as shown in Figs. 
4 and 5, the trapped ion cannot be a good prober and a cooler, because the phonon exchange is not 
complete between the ion and the cantilever. Both calculations with and without the RWA clearly 
demonstrate this point, while the non-RWA treatment could give a more accurate description of 
quantum behavior of the system. 

IV Discussion and conclusion 

One point we should mention is that, we have only studied the first step in the proposal [14], 
because the second and the third steps have been well investigated by standard approaches. 
According to [14], after the vibrational quanta has been fully transferred to the ion, we switch off 
the coupling K , and meanwhile turn on the red- or blue -detuned laser radiation. The current ion 
trap techniques have successfully coupled the internal and the vibrational degrees of freedom of 
the trapped ultracold ion [18]. But we have not yet found any experimental report for detecting the 
vibrational quanta of an ion beyond the sideband cooling regime. Considering an almost full 

phonon transfer from the cantilever to the ion, we choose n a (0) to be 6 so that the vibrational 

quanta obtained by the ion can be detected at the level of a single quanta by means of the 
sophisticated optic method [15,16]. Nevertheless, from [14], it seems to be not difficult 

experimentally to detect a cantilever with the vibrational quanta n a (0) > 4000 . As it in 

principle works for any n a and n b , our method, providing an exact solution for the evolution 
of system, could be very helpful for the future experiment. 

On the other hand, however, the phenomenological treatment of the decay in Eq. (1) due to 
thermal environment restricts our approach to work only on the cases of few vibrational quanta. 
Strictly speaking, if the vibrational quanta owned by the cantilever and by the trapped ion are not 
approaching zero, we cannot simply write the vibrational damping to be like in Eq. (1). Instead, 
we should utilize density matrix, not simply wave functions, to study the dynamics of the system 
[19]. Nevertheless, under the simple model with an ultracold ion coupled to a cantilever with only 
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6 vibrational quanta, Eq. (1) is valid in the case of max{r a , T b } « K . This approximation goes 
better when /f/maxjr^r^Jis larger. In this sense, our investigation of the quantum dynamics 

of the system is exact only within the considered simple model with few vibrational quanta 
involved. 

Besides, we have to mention again that Eq. (1) itself is an approximate description of the real 
coupling system without including three - and higher - order perturbative expansion terms of the 
Coulomb energy. So the model under consideration only works for the weak coupling cases. If we 
want to investigate a very strong coupling model, Eq. (1) should be modified to involve 
multi-phonon terms due to the higher-order terms of the perturbative expansion, which would be 
an interesting work in the future. What we have done here is to present an exact quantum behavior 
of the coupling system within the simple model in cases of relatively strong couplings (i.e., 1.9 
MHz < K < 5.0 MHz in our solution). 

To show the non-RWA effect experimentally, we need large coupling constant K . As defined in 

[14], JCocV Q /(mMVG)d 6 ) 1 ' 2 where m and M are the mass of the ion and the cantilever, 

respectively, and Vo and d are defined in Fig. 1. In terms of [3,14], the vibrational frequencies of 
the cantilever and the ion can be 19.7 MHz, and the available coupling constant K is about 0.3 
MHz. So if we decrease the vibrational frequencies of both the cantilever and the ion by, for 
example, 4 times, K would be about 1.2 MHz, weaker than other two characteristic frequencies 
by only 4 times. As a result, the non-RWA effect studied in this paper would appear. Alternatively, 
the choice of a lighter cantilever, e.g., a single wall carbon nanotube, could also increase the 
coupling K . 

In summary, we have studied the time evolution of an electrostatically coupling system involving 
a trapped ultracold ion and a precooled cantilever, which is the first step of a previous proposal to 
employ the trapped ion as a transducer and to fast cool the cantilever down to its ground 
vibrational state. With the squeezed-state expansion technique, we can exactly demonstrate the 
quantum behavior of two coupled subsystems within the model given in [14], which is very 
different from the solution with the RWA under certain conditions. Since the field of QEMS is 
developing very quickly and is playing more and more important role in QIP for the scalability [10] 
and the qubit readout [11-13], we believe that our method as well as the results presented here 
could be useful in the fields of QEMS and QIR 
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Appendix I 

From the initial conditions CC 2 (0) = 0, CC 3 (0) = , and ^£=0|f=0) = l,we have 

|p(0)| 2 (0|e a ' (0)aa e a ' (f)a+a+ |0) 

= \p(0)\ 2 {0\Y J -[a l (0)a 2 T ■ E^i(0)(^ + ) 2 r|0> 
n' rrt 

is— i f j 

n * m 

= |p(o)| 2 (oe^[«(o))^ 2/i ] ■ ^mm'iayp) 

= \P(°f (° E 777" (0)| 2 "(2«) ! ]| 0). 
'TOO 

It is easy to prove 

^1-41^(0)1 

So we get . — l Similarly, from (t = |«3 + «2| f = 0) = iV a , we can obtain 

A /(l-4|^ 1 (0)| 2 ) 

^ (0)2 -i = K, 



(l-4|a,(0)| 2 )^ 



where we have used 



n=o(n\) (1-4^(0) \Y 2 



Appendix II 

For solving Eqs. (13) - (15), we set a x (f) = a w e at , a 2 (t) = a 2Q e at and a 3 (t) = a 3Q e at , 
and obtain 

£la w = 2{to-iT a )a w - koc 2Q 

j£la 20 =(co+v - iT a - iT b )a 2Q - 2jca lQ - 2jca 30 

£la 3Q =2(v-iT h )a 30 -Ka 20 . 

By diagolizing the related determinant, we get to 
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a: 



,(1,2,3) 
'20 



2{(0-iT a )-£l 



a 



,(1,2,3) 
10 



K 



a 



,(1,2,3) 
'30 



2(G)-ir b )-jn (W) 



or. 



.(1,2,3) 
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where 



a (I) =-[r a +r b +j(a)+v),] 

n (2) = -vr a + r b + j(o)+v)] + V-[y(r a -rj-«+v] 2 -r 2 , 
a (3) = -[n, + r, + i{co+ v)] - ^-u{T a -T b )-o)+vf -k 2 . 

By means of the initial conditions 0C X (0) = ^ ->[Nj{N a + 1), ^ 2 (0) = 0> a 3 (0) = , we obtain 
Eqs. (16) and (17). 



Fig. 1 Schemetic diagram of a single trapped ultracold ion coupled to a doubly clamped cantilever, 
where d is the separation comparable to the vibrational deviations of the ion and the cantilever 
from their equilibrium positions, and V is the bias voltage. 



Fig.2 Time evolution of n a and n b calculated under and beyond the RWA, where co=19.7 MHz , 

v=19.7 MHz , k=1.8 MHz , I>r b =0.0197 MHz, and 77 a (0) = 6,^(0) = . 

Fig. 3 Time evolution of n a and n b calculated under and beyond the RWA, where co=19.7 MHz , 

v=19.7 MHz , k=5.0 MHz , I>r b =0.0197 MHz, and n a (0) = 6,n b (0) = . 

Fig.4 Time evolution of n a and n b calculated under and beyond the RWA, where co=19.7 MHz , 

v=16.0 MHz , k= 4.0 MHz, I>r b =0.0197MHz, and n a (0) = 6,n b (0) = . 

Fig.5 Time evolution of n a and n b calculated under and beyond the RWA, where co=19.7 MHz, 

v=16.0 MHz , k=5.0 MHz, r a =r b =0.0197 MHz, and n a (0) = 6,n b (0) = . 



Captions of the Figures 
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